We characterize the importance of metallicity on the pres- 
ence of molecular hydrogen in damped Lyman-a (DLA) sys- 
tems. We construct a representative sample of 18 DLA/sub- 
DLA systems with log N(U i) > 19.5 at high redshift (z abs > 1 .8) 
with metallicities relative to solar [X/H] > -1.3 (with [X/H] = 
log JV(X)/jV(H)-log(X/H) and X either Zn, S or Si). We gather 
data covering the expected wavelength range of redshifted H2 
absorption lines on all systems in the sample from either the 
literature (10 DLAs), the UVES-archive or new VLT-UVES 
observations for four of them. The sample is large enough to 
discuss for the first time the importance of metallicity as a cri- 
terion for the presence of molecular hydrogen in the neutral 
phase at high-z. From the new observations, we report two new 
detections of molecular hydrogen in the systems at z a bs = 2.43 1 

^ toward Q 2343+125 and z abs = 2.426 toward Q 2348-011. 

f--*) ■ We compare the H2 detection fraction in the high-metallicity 

£vq sample with the detection fraction in the overall sample from 
1 Ledoux et al. (2003). We show that the fraction of DLA systems 
£3 i with log / = log 2AT(H 2 )/(2jV(H 2 ) + N(U 1)) > -4 is as large as 

^ 50% for [X/H] > -0.7 when it is only -5% for [X/H] < -1.3 
and -15% in the overall sample (with -2.5 < [X/H] < -0.3). 

,— ' This demonstrates that the presence of molecular hydrogen at 
high redshift is strongly correlated with metallicity. 

> ■ 

^sO [ Key words, galaxies: ISM - quasars: absorption lines - 
t> ■ quasars: individual: Q 2343+125; Q 2348-011 
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Abstract. 



1. Introduction 

Early searches for molecular hydrogen in DLAs, though not 
systematic, have lead to either small values of or upper limits 
on the molecular fraction of the gas. For a long time, only the 
DLA at z a b s = 2.811 toward Q 0528-250 was known to con- 
tain H2 molecules (Levshakov & Varshalovich 1985, Foltz et 
al. 1988). Ge & Bechtold (1999) searched for H 2 in a sample 
of eight DLAs using the MMT moderate-resolution spectro- 
graph (FWHM = 1 A). Apart from the detection of molecu- 
lar hydrogen at z- d bs = 1-973 and 2.338 toward, respectively, 
Q 0013-004 and Q 1232+082, they measured in the other sys- 
tems upper limits on / in the range 10~ 6 - 10~ 4 . 

A major step forward in understanding the nature of 
DLAs through their molecular hydrogen content has recently 
been made possible by the unique high-resolution and blue- 
sensitivity capabilities of UVES at the VLT. In the course of 
the first large and systematic survey for H2 at high redshift, we 
have searched for H2 in DLAs down to a detection limit of typ- 
ically AT(H 2 ) = 2 x 10 14 cm" 2 (Petitjean et al. 2000, Ledoux 
et al. 2003). Out of the 33 surveyed systems, eight had firm 
detections of associated H 2 absorption lines. Considering that 
three detections were already known from past searches, H 2 
was detected in ~ 15% of the surveyed systems. The existence 
of a correlation between metallicity and depletion factor, mea- 
sured as [X/Fe] (with X = Zn, S or Si) was demonstrated (see 
also Ledoux et al. 2002a) and the DLA and sub-DLA systems 
where H 2 was detected were usually among those having the 
highest metallicities. 



* Based on observations carried out at the European Southern 
Observatory (ESO) under prog. ID No. 65.O-0158, 65.O-0411, 67.A- 
0022, 67.A-0146, 69.A-0204, 70.B-0258, 072.A-0346 with UVES in- 
stalled at the Very Large Telescope (VLT) on Cerro Paranal, Chile 



However, the high metallicity end of the Ledoux's sam- 
ple was biased by the presence of already known detections. 
Therefore, to investigate further this possible dependence with 
metallicity, and derive what is the actual molecular content of 
the high-redshift gas with highest metallicity, we have searched 
a representative sample of high metallicity DLAs for H 2 . The 
number of H 2 measurements in systems with [X/H] > -1.3 is 
twice larger in our sample compared to previous surveys. We 
describe the sample and the observations in Section 2, present 
two new detections of H 2 in Section 3 and the results of the 
survey and our conclusions in Section 4. 

2. Sample and Observations 

We have selected from the literature all z > 1.8 DLAs and sub- 
DLAs systems (logA^(Hi) > 19.5) with previously measured 
elemental abundances (e.g. Prochaska & Wolfe 2001, Kulkarni 
&Fall2002) larger than [X/H] > -1.3 and accessible to UVES. 
The inclusion of sub-DLAs is justified by the fact that for 
logV(Hi) > 19.5 most of the hydrogen is neutral (Viegas 
1995). In addition, Ledoux et al. (2003) have shown that for 
log N(H 1) > 19.5, there is no correlation between the presence 
of H 2 and the H 1 column density. The fact that we include ALL 
known systems with these criteria guarantees that the sample is 
representative of the population of DLA-subDLA. 

We ended up with a sample of 15 DLAs and 3 sub- 
DLAs (with log N(H 1) = 19.7, 20.10 and 20.25) pre- 
sented in Table 1. The H 2 content of ten of these sys- 
tems had already been published (Srianand & Petitjean 2001, 
Petitjean et al. 2002, Ledoux et al. 2002b, 2003, 2006b, 
Srianand et al. 2005 and Heinmuller et al. 2006). Five of 
the eight remaining systems had data in the UVES archive 
(Q 1209+093: Prog. 67.A-0146 P.I. Vladilo; Q 21 16-358: 
Prog. 65.O-0158 P.I. Pettini; Q 2230+025: Prog. 70.B-0258, 
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P.I. Dessauges-Zavadsky; Q 2243-605: Prog. 65.O-0411 P.I. 
Lopez; Q 2343+125: Prog. 69.A-0204 and 67.A-0022, P.I. 
D'Odorico). New observations of Q 0216+080, Q 2206-199, 
Q 2343+ 125 and Q 2348-011 have been performed with the 
Ultraviolet and Visible Echelle Spectrograph (UVES, Dekker 
et al. 2000) mounted on the ESO Kueyen VLT-UT2 8.2 m tele- 
scope on Cerro Paranal in Chile. These new observations re- 
sulted in two new detections described in the next Section. 
The data for each of the eight QSOs were reduced using the 
UVES pipeline which is available as a context of the ESO 
MIDAS data reduction system (see e.g. Ledoux et al. 2003 for 
details). Standard Voigt-profile fitting methods were used for 
the analysis of metal lines and molecular lines, when detected, 
to determine column densities using the oscillator strengths 
compiled in Ledoux et al. (2003) for metal species and the os- 
cillator strengths given by Morton & Dinerstein (1976) for H2. 
We adopted the Solar abundances from Morton (2003) based 
on meteoritic data from Grevesse & Sauval (2002). 

The characteristics of the sample are summarized in 
Table [0 The H 1 column densities and most of metallicities 
are from the compilation by Ledoux et al. (2006a). Slight dif- 
ferences with previously published values, e.g. Ledoux et al. 
(2003), are due to the use of different damping coefficients for 
H 1. The only known z > 1.8 H2 bearing DLA system out of 
this sample is the system at z a b s = 2.337 toward Q 1232+082 
(Srianand et al. 2000). 

3. Two new detections 

We report two new detections of H2 in DLAs from our new 
observations. Detail analysis and interpretation of the physical 
conditions in these two DLAs are out of the scope of the present 
work and will be described in a subsequent paper. 
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Fig. 1. Velocity plots of observed H2 absorption lines from 
the J = and J = 1 rotational levels at z a t, s = 2.43127 to- 
ward Q 2343+125. Profiles of associated Si 11 and Fe 11 ab- 
sorptions are shown in the upper panels. The best-fitting model 
to the system is overplotted; the location of Voigt-profile sub- 
components are shown as vertical dashed lines. 



3.1. Q2343+125,z abs = 2A3\ 

This DLA system has first been studied by Sargent et al. 
(1988). High resolution data have been described by Lu et 
al. (1996), D'Odorico et al. (2002) and Dessauges-Zavadsky 
et al. (2004). The profile of the metal lines is spread over 
more than 250 km s" 1 from z abs = 2.4283 to 2.4313 but the 
strongest component is centered at z a b s = 2.43127 correspond- 
ing to the red edge of the above redshift range. From Voigt- 
profile fitting to the H 1 Lyman-a, j3 and y lines, we find that 
the damped Lyman-a line is centered at z a b s = 2.431 and the 
column density is log N(H 1) = 20.40 ± 0.07, consistent with 
previous measurement by D'Odorico et al. (2002; log N(H 1) = 
20.35 ± 0.05). We use Zn as the reference species for metallic- 
ity measurement and find [Zn/H] = -0.89+0.08. This is con- 
sistent with previous findings. Absorption from the J = 1 and 
probably from the J = rotational levels of H2 is detected in 
this system at z a b s = 2.43127 (see Fig.^. The optically thin 
H2 absorption lines are very weak, i.e. close to but above the 
3cr detection limit. A very careful normalization of the spec- 
trum has been performed, adjusting the continuum while fitting 
the lines. The best-fitting consistent model for H2 is shown in 
Fig-H The total H2 column density, integrated over the J = 
and 1 levels, is estimated to be log N(H 2 ) = 13.69 ± 0.09 



(12.97+0.04 and 13. 60+0. 10 for J = and 1 respectively). This 
leads to the smallest molecular fraction observed up to now, 
log / = -6.41^j [g. We also derive an upper limit on the detec- 
tion of absorption from the J = 2 level, log Af(H2-J=2) < 13.1 
at the 3<r level. 

3.2. Q 2348-0 11, = 2.426 

There are two DLA systems at z a b s = 2.426 and z a b s = 2.615 
toward Q 2348-011, with total neutral hydrogen column den- 
sities of respectively log N(U 1) = 20.50 ±0.10 and log N(U 1) 
= 21.30 ± 0.08. Conspicuous H2 absorptions are detected in 
the z a bs — 2.426 DLA system, the only system to be consid- 
ered here as displaying a high metallicity (see Fig. |3. The 
molecular lines are very numerous and strong but the spec- 
tral resolution of our data is high enough to allow unambigu- 
ous detection and accurate determination of the line parame- 
ters. Seven H2 components spread over about 300 km s _1 were 
used for the H2 fit. It is interesting to note that the strongest 
metal component (at V = km s _I in Fig. [3 has no associated 
H2 absorption. Strong H2 absorption is seen at V ~ -150 and 
+50 km s . All seven molecular components have associated 
C 1 absorption. However, additional components are needed to 
fit the metal absorption lines: 9 components for C 1 and 13 com- 
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Table 1. Metal and molecular content of high-metallicity DLA/sub-DLA systems at tabs > 1.8 
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1 Neutral hydrogen column densities and metallicities relative to Solar, [X/H] = log[7V(X)/W(H)] 
reference element when Zn 11 is detected, or else either S or Si) are from Ledoux et al. (2006a). 



log[/V(X)/iV(H)] e (with X = Zn as the 

and J = 1 are given in 



2 Molecular hydrogen column densities are summed up over all / levels in case of detection and upper limits for J 
case of non-detection. 

3 Molecular fraction, / = 22V(H 2 )/(2iV(H 2 ) + A'(Hi)). 

a Petitjean et al. (2002), b Ledoux et al. (2003), c This work, d Srianand et al. (2005), e Heinmiiller et al. (2006), / Ledoux et al. (2002b), g 
Srianand & Petitjean (2001), h Ledoux et al. (2006b) 
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Fig. 2. Absorption profiles at z a b s = 2.4263 (taken as zero of the 
velocity scale) toward Q 2348-01 1 . The positions of the seven 
(respectively 13) components needed to fit the H2 (respectively 
metal line) profiles are indicated by vertical dashed lines. The 
resulting best-fitting model to the system is overplotted. 



ponents for the singly ionized species. H2 absorption from the 
rotational levels J = to 5 are unambiguously detected. The 
total H2 column density integrated over all rotational levels is 
log ./V(H2) = 18.45^ 2g, corresponding to a molecular fraction 
log / = - 1 .76^Q3g. We also derive an upper limit on the column 
density of HD molecules, leading to log A^(HD)/A^(H 2 ) < -3.3. 



4. Metallicity as a criterion for the presence of 
molecular hydrogen 

From Table 1, it can be seen that H2 is detected in nine high- 
metallicity systems out of 18. However, for two of the detec- 
tions, the corresponding value of log / is smaller than most of 
the upper limits derived for other systems. All upper limits are 
smaller than -4.5 and all detections larger than these upper lim- 
its are larger than -3. We therefore define a system with large 
(respectively small) H2 content if log / is larger (respectively 
smaller) than -4. 

In Fig. we plot the molecular fraction, log /, versus 
metallicity, [X/H], for our representative sample of DLAs with 
[X/H] > -1.3 (18 measurements summarized in Table 1) and 
measurements by Ledoux et al. (2003) for [X/H] < -1.3 (23 
measurements). The log / distribution is bimodal with an ap- 
parent gap in the range -5 < log / < -3.5 justifying the above 
classification of systems. Note that this jump in log / has al- 
ready been noticed before by Ledoux et al. (2003) and is sim- 
ilar to what is seen in our Galaxy (Savage et al. 1977; see also 
Srianand et al. 2005). It is apparent that the fraction of sys- 
tems with molecular fraction log / > -4 increases with in- 
creasing metallicity. It is only ~5% for [X/H] < -1.3 when it is 
-39% for [X/H] > -1.3. This fraction is even larger, 50%, for 
[X/H] > -0.7 which is the median metallicity for systems with 
[X/H] > -1.3. In addition, all systems with [X/H] < -1.5 have 
log /< -4.5. 

We conclude that metallicity is an important criterion for 
the presence of molecular hydrogen in DLAs. This may not 



P. Petitjean et al.: H2 in DLA systems 



5 



1 1 1 1 1 , 1 1 1 1 , 

-f- 


1 1 , 


1 , 1 1 
1 — 1 1 — 1 


\ .... \ .... \ . 


bD— 1 1 iTi 
, , 1 . , 


- 
- 

. 1 . , , , 



I I I I I I I I I I I I I I I I I I I I I I I I I I 

-2.5 -2.0 -1.5 -1.0 -0.5 0.0 

[X/H] 



Fig. 3. Logarithm of the molecular fraction, 
/ = 2N(H 2 )/(2N(H 2 ) + N(H 1)), versus metallicities, [X/H] = 
log 7V(X)//V(H)-log(X/H) and X either Zn, S or Si, in DLAs 
from the sample described in this paper ([X/H] > -1.3, see 
Table 1) and the sample of Ledoux et al. (2003, [X/H] < -1.3). 
Filled squares indicate systems in which H2 is detected. 
Dashed lines indicate the limits used in the text (log / = -4; 
[X/H] = —1.3). The dotted line indicate the median of the 
high-metallicity sample ([X/H] = -0.7). 



be surprizing as the correlation between metallicity and deple- 
tion of metals onto dust grains (Ledoux et al. 2003) implies 
that larger metallicity means larger dust content and therefore 
larger H2 formation rate. In addition, the presence of dust im- 
plies a larger absorption of UV photons that usually dissoci- 
ate the molecule. More generally, Ledoux et al. (2006a) have 
shown that a correlation exists between metallicity and veloc- 
ity width in DLAs. If the latter kinematic parameter is inter- 
preted as reflecting the mass of the DM halo associated with 
the absorbing object, then DLAs with higher metallicity are as- 
sociated with objects of larger mass in which star formation 
could be enhanced. All this makes it arguable that, in DLAs, 
star-formation activity is probably correlated with the molecu- 
lar fraction (Hirashita & Ferrara 2005). It is therefore of first 
importance to survey a large number of DLA systems to de- 
fine better their molecular content and use this information to 
derive the physical properties of the gas and the amount of star- 
formation occuring in the associated objects. 
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